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Abstract 

The equilibrium thermodynamic properties of a linear polymer chain confined to a space between 
two impenetrable walls (lines) at a distance D under various solvent conditions have been studied 
using series analysis and exact enumeration technique. We have calculated the end to end distance 



of polymer chain, which shows a non-monotonic behaviour with inter wall separation D. The 
density distribution profile shows a maxima at a particular value of {D =)D* . Around this D* , 
our results show that the collapse transition occurs at higher temperature as compared to its bulk 

value of 2d and 3d. The variation of 9— temperature with D shows a re-entrance behaviour. We 

< i 

iL? ■ also calculate the force of compression exerted by the walls (lines) on the polymer. 
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I. INTRODUCTION 

In last few years 1-6 much attention has been paid on the subject of confinement of polymer 
and its scaling relations. This is due to the fact that a polymer chain in restricted geometries 
exhibits specific and interesting properties which find application in steric stabilization of 
colloidal dispersions, thin films, adsorption behaviour of gels, surface coatings, membrane in 
nano pores etc 6 ' 7 . Several theoretical and numerical attempts has been made to understand 
these properties 1 ' 6,10 A polymer chain is said to be confined when the spacing D between the 
parallel walls (lines in 2d) is less than the average end to end distance < Re > of the chain. 
The parallel (< Re\\ >)and perpendicular component (< Re± >) of < R E > is expected to 
behave as 2 , in 3d, 

3D 3 

< R E \\ >~a.A^4(— ) 1 "4i7 (1) 

(Jj 

and 

< R E ± >~ D. (2) 

Where N is the total number of monomers in the chain, a is the lattice parameter and v is 
the end to end distance exponent respectively. When D is greater than < Re >, the effect 
of confinement on its conformational properties vanishes i.e. 

< R E{{ > = < R E± >~ a N u (3) 

The monomer density (p(z)) profile near the wall is expected to behave as 

p(z) ~ Z" (4) 



here z is the perpendicular distance from the surface. It is assumed that the force exerted 
by the polymer on the wall is proportional to the monomer density near the wall. It has 
been shown that 9 

hm_o^ ~ B-L- (5) 

where k is the non-universal amplitude which relates the end to end distance < Re > of a 
polymer to the chain length N, B is the universal amplitude ratio, / is the force exerted on 
the wall due to confinement of the chain and T is the temperature of the system respectively. 



Considerable efforts have been made 1-10 to verify the scaling relations and determine the 
value of B in two and three dimensions. In a recent simulation by Hsu and Grassberger 1 
much better value of 5(=2.04 ± 0.04 (in 2d) and 2.13 ± 0.11 (in 3d)) has been obtained as 
compared to the previous Monte Carlo estimates 2 ' 10 . 

The complete phase diagram of surface interacting polymer chain in poor solvent con- 
dition has been studied recently 11 ' 12 . Our aim in this paper is to analyze the effect of 
geometrical confinement such as the one given by a slit (or strip) on coil-globule transition 
under good and poor solvent conditions. In any Monte Carlo simulation, and particularly 
in the study of polymeric system, improper sampling of hyper-space under finite number of 
steps being considered which some time overlook some of the interesting features of these 
systems 11 ' 13 . Keeping this in mind, we have analyzed the partition function through series 
analysis and exact enumeration technique. We prefer this technique because in this case the 
complete partition function can be analyzed exactly and the scaling corrections are correctly 
taken into account by a suitable extrapolation scheme 13 ' 14 . For the sake of simplicity and 
comparison with the previous results the walls here have been considered as neutral. 

The work is organized as follows; In Section 2, we briefly describe the model and the 
method. Section 3 deals with the monomer density profile in two and three dimensions 
where we compare our results with the known results. In Section 4 we study the effect of 
confinement on the coil-globule transition and establish the phase diagram which shows the 
reentrance behaviour in three dimensions. The paper ends with the discussion in Section 5. 

II. MODEL AND METHOD 

We represent a polymer chain by self-avoiding walks (SAWs) on a square (2d) and cubic 
(3d) lattices. The polymer chain is confined in between strip of two impenetrable surface of 
height D = 0,1, 2. .16 in case of square lattice and D = 0,1, 2. .8 for cubic lattice, respectively. 
The walk originates from the middle of the surface (kept fixed) and spreads along all possible 
directions. This ensures that the chain is not affected by the presence of the upper surface 
for large values of D and one recovers the usual statistics of SAWs. The thermodynamic 
properties associated with the confined polymer may be obtained from the partition function 



which can be written as a sum over all possible configurations confined in distance D i.e. 



Z N {D) = Y,C N {D) 



X 



N 



(6) 



here x is the fugacity associated with each step of the walk and Cn{D) is the number distinct 
SAWs of N steps confined in width D. We have obtained Cn{D) up to N < 28 in two 
dimensions (for D < 16) and N < 19 (for D < 8) in three dimensions. In general, it is 
appropriate to assume that for large N, 



C N {D) ~ AP"-V(£) 



iY 



(7) 



here n(D) is the connectivity constant of the lattice and 7 is the associated critical exponent. 
The value of n{D) can be estimated using the ratio method 13-14 with an associated Neville 
table. These values are found to be in good agreement with the result obtained by the exact 
transfer matrix calculation 15 . The end to end distance exponent v may be found using the 
relation 

log(< R N > J < R N _ 2 >) 



v 



(8) 



log(iV/iV - 2) 

where < Rn > is the end to end distance exponent of iV steps walk. The values of v — 
0.748 ± 0.002(2<i) and 0.58 ± 0.01(3<i) are in good agreement with the known results in 
the limit D — > oo 6 . For finite value of D, we restore the scaling proposed by Milchev and 
Binder 2 . The slopes found for 2d and 3d for parallel and perpendicular component of end to 
end distance are in good agreement with Eqs. (1 — 2) as shown in Figs. (1 — 2). However, 
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FIG. 1: Scaling for parallel component of end to end distance exponent. 
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FIG. 2: Scaling for perpendicular component of end to end distance. 

the variation < Re > is non-monotonic with D as shown in Fig. 3. We find that < Re > 
first decreases with decreasing D before it starts to rise. This behaviour was ascribed earlier 
to a lateral squeezing of chains with increasing D, followed by an elongation in the direction 
parallel to the walls 16-17 
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FIG. 3: The variation of end to end distance (< Re >) of polymer chain for various length with 
inter wall separation D for 3d and 2d. 



III. MONOMER DENSITY DISTRIBUTIONS 



In the present model, one end of the chain is anchored at the middle of one of the 
walls. Consequently, it introduces asymmetry in the distribution of monomer density in 
the direction perpendicular to the plane of wall. In earlier simulations 1 , the walker has 



been allowed to start from any of the planes in between the walls, therefore, the monomer 
distribution in such case is found to be symmetric. The density distribution of polymer 
attached to the surface has been shown in Fig. 4 for 2d and 3d. From these figures it is 
obvious that when the thickness of the strip decreases, the monomer density near the wall 
increases. In Fig. 5, we show the monomer density for fixed D which increases with iV,and 
therefore, is not a finite size effect. In fact the peak shifts away from surface, at which 
polymer is grafted, as N increases. 
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FIG. 4: The variation of monomer density with D for different value of D. 

The force exerted on the wall is proportional to the monomer density on the surface and 
may be expressed in terms of the work done by moving one of the surfaces 
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Using the simple scaling ansatz 



dG k B T dfj, 
dD ~ fi dD' 



P(z) 
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-^(i-0 



l/Vd 



(9) 



(10) 



where ^ = y/D + 1 and v& is the end to end distance exponent of d dimension, with 
A = 10.38(2d) and 18.74(3d), as proposed by Hsu and Grassberger 1 we find the value of the 
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universal amplitude ratio B to be equal to 2.01(2d) and 2.1 (3d), respectively. These values 
are in good agreement with the exact result 9 and simulations 1 . 
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FIG. 5: The density profile with D for different chain length. These figures show that the peak 
height increases with N but shifts out to larger D as N increases. 



IV. EFFECT OF CONFINEMENT ON COIL-GLOBULE PHASE TRANSITION 



As confinement affects the monomer density profile, one expects a significant change in 
the transition point if one varies the distance between the two walls. A square lattice with 
finite D is essentially an one dimensional system and hence there is no phase transition as 
such. Keeping this in mind, we study such an effect, on a simple cubic lattice with variable 
wall thickness D. We consider self- attracting-self avoiding walks (SASAWs) and associate 
energy — e p for nearest neighbour monomers which are not chemically bonded. The partition 
function of the confined chain may be written as 

[C N , Np (D)u N * (11) 



Z N (D,u) 



N v 



where u is the Boltzmann weight defined as u = e e p/ k P T and N p is the number of pairs 
of monomer of chain length N. Cn,n p (D) is the number of distinct configurations with N p 



number of nearest neighbour pairs confined between the walls at distance D. The reduced 
free energy can be written as 
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FIG. 6: The variation of u with D shows re-entrance behaviour in 3d. The dip in u indicates the 
rise in 6— temperature. It may be noted that this minimum value of u c shown in phase diagram 
is less than bulk (u c =1.76 in 3d) and surface (« c =1.93 in 2d) value. 



G(D,u) = lim — log Z N (D,u) = log[/i(L>, «)]• 



(12) 



To characterize conformational change in the system, we study specific heat defined as 
C v = d 2 G/dei as a function of temperature. The peak in the C v — T diagram has been 
identified as a sudden change in conformation. Such changes in heat capacity in a narrow 
temperature range are due to large fluctuations in energy. Studying the distribution of end to 
end distance in this region we identify such transition as Coil -Globule transition, ofcourse, 
a sharp transition of a single chain occurs in the limit N — > oo only, and hence the accuracy 
of the extrapolation in Eq. 12 is crucial. In the limit D — > oo, we find u c = 1.93 and 1.76 for 
2d and 3d, respectively. These values are in good agreement with earlier results 18 " 20 . The 
variation of u c , with D has been shown in Fig. (6) which shows a re-entrance behaviour. 
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It can be seen from Fig. (6) that at a fixed value of monomer-monomer attraction with a 
change of inter wall separation D, a polymer which was initially in the coil phase goes into 
the globule phase. With further increase of D the polymer again attains the coil phase. 

V. CONCLUSIONS 

In this paper, we have studied the effect of confinement on the equilibrium thermodynamic 
properties of a linear polymer chain in various solvent conditions. The scaling relations found 
by us for < Ru > and < R± > are in agreement with earlier predictions 2 . We show that due 
to the confinement, the end-to-end distance, < ReI varies non-monotonically with inter- 
wall separation D. This is also reflected in variation of monomer density with maxima at 
particular D*. In this region the number of non-bonded nearest neighbours increases due 
to confinement and hence there is a rise in the transition temperature (i.e. decrease in u c ) 
as shown in Fig. (6). The re-entrant behaviour can also be seen from Fig. (6). It has 
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FIG. 7: The variation of force exerted by confining walls on per monomer of the polymer chain 
with D. 



been observed that for a given value of monomer-monomer attraction a polymer which was 
initially in the coil phase goes into the globule phase due to the rise in nearest neighbour 



pairs due to confinement at the particular D* . With further increase or decrease of the 
wall's separation D, the number of these pairs decreases and hence the polymer acquires the 
coil phase. 

We also studied this effect on a square lattice with finite D. We do find a peak in the 
specific heat which signals that the system has a tendency to be in the ordered state, but this 
tendency will only turn into a phase transition when the slit width diverges i.e. D > N u . 

It is interesting to see the variation of force exerted by the confining wall with D in Fig. 
(7). It is qualitatively similar to the one obtained in a simulation 10 , in which the force has 
been applied by atomic force microscope tip on polymer chain attached to the surface. In 
such a case it has been found that polymer escapes from an approaching AFM tip (finite 
surface)by forming a stretched tether 10 . Since in our work the surface has been considered 
as an infinite plane, "escape transition" is not possible. However, the work with a finite 
surface of specific geometry is in progress. 
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